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Graphene is famous for being a host of 2D Dirac fermions. However, spin-orbit coupling intro¬ 
duces a small gap, so that graphene is formally a quantum spin Hall insulator. Here we present 
symmetry-protected 2D Dirac semimetals, which feature Dirac cones at high-symmetry points that 
are not gapped by spin-orbit interactions and exhibit behavior distinct from both graphene and 
3D Dirac semimetals. Using a two-site tight-binding model, we construct representatives of three 
possible distinct Dirac semimetal phases, and show that single symmetry-protected Dirac points are 
impossible in two dimensions. An essential role is played by the presence of nonsymmorphic space 
group symmetries. We argue that these symmetries tune the system to the boundary between a 
2D topological and trivial insulator. By breaking the symmetries we are able to access trivial and 
topological insulators as well as Weyl semimetal phases. 


Over the past decade, graphene has attracted intense 
interest as a material with Dirac cones at the Fermi en¬ 
ergy and, as a consequence, a number of unique electronic 
properties [1, 2]. The Dirac points in graphene, as in sim¬ 
ilar materials [3-6], are protected by symmetry, but only 
in the absence of spin-orbit coupling. Spin-orbit coupling 
opens a gap at the Dirac point, leading to a topologi¬ 
cal insulating phase [7, 8]. The discovery of topological 
insulators (TIs) heightened interest in three-dimensional 
Dirac semimetals, which host 3D Dirac points when spin- 
orbit coupling is included [9, 10]. The concept of Dirac 
and Weyl superconductors has also been recently intro¬ 
duced [11]. In this Letter we introduce a system that has 
symmetry-protected 2D Dirac points in the presence of 
spin-orbit coupling and provide a classification of such 
systems in general. These are of interest because they 
are symmetry tuned to the boundary between topologi¬ 
cal and trivial insulating phases. 

Three-dimensional Dirac semimetals fall into two dis¬ 
tinct classes. In Ref. [10] we introduced a Dirac 
semimetal with Dirac points at high-symmetry points on 
the surface of the Brillouin zone (BZ). Candidate mate¬ 
rials include /?-cristobalite Bi02 [10], as well as distorted 
spinel materials such as BiZnSi04 [12]. In these materials 
the semimetallic state is at the boundary between strong 
and weak topological insulating phases, and an essential 
role is played by the nonsymmorphic symmetry of the 
crystal space group. A distinct class of Dirac semimetals 
was introduced in Refs. [13, 14], and has been observed in 
Cd3As2 and NasBi [15-17]. Here, the Dirac points arise 
due to a band inversion and occur at a generic point 
on a C3 symmetry axis in the interior of the BZ. Open¬ 
ing a gap by lowering the symmetry in these materials 
necessarily leads to a topological insulator - the trivial 
insulator is not adjacent. The 2D Dirac semimetals we 
introduce here are analogous to the former class: they 
arise due to a nonsymmorphic symmetry that requires 
the conduction and valence bands to touch and exist in 


the presence of significant spin-orbit coupling. We will 
argue that the nonsymmorphic character correlates with 
the fact that they are at the boundary between a triv¬ 
ial and topological insulator. We will begin by clarify¬ 
ing the role of nonsymmorphic symmetries in protecting 
Dirac points. We will then introduce a simple model 
system that allows us to characterize the allowed Dirac 
phases in 2D, and conclude with a brief discussion of the 
possible material venues for these phases, including the 
layered iridium oxide superlattices recently proposed and 
studied in Ref. [18]. 

It has long been known that nonsymmorphic symme¬ 
tries lead to extra degeneracies in electronic band struc¬ 
tures that cause bands to “stick together” due to the ex¬ 
istence of higher-dimensional projective representations 
of the little groups of certain values of k[I9]. This fact 
can be understood as a simple consequence of fractional 
translation symmetries. Nonsymmorphic space groups 
are distinguished by the existence of symmetry opera¬ 
tions that combine point group operations g with trans¬ 
lations t that are a fraction of a Bravais lattice vector. 
In 2D, the relevant operations, denoted {g|t}, are screw 
axes g = C'2ri^ in± T z), glide mirror lines g = Mn^ and 
glide mirror planes g = M^, in conjunction with a half¬ 
translation t that satisfies gt = t as well as * = — I 
for the “odd” reciprocal lattice vectors G. The con¬ 
sequence of the fractional translation symmetry is sim¬ 
plest in the case where the unit cell of simple symmor- 
phic crystal is doubled. In that case, the folded back 
bands necessarily touch on lines in momentum space. 
If the symmetry is reduced so that {A|t} is violated, 
the degeneracy is, in general, split. However a non¬ 
symmorphic symmetry {git} still protects degeneracies 
in the invariant line or plane in the BZ that satisfies 
gk = k. In this invariant space the Bloch states can be 
chosen to be eigenstates {g|t}|M^) = ±Ae**'*|u^). For 
k — >• k -I- G, with * = — I, the two eigenstates must 
switch places [20]. In the absence of other degeneracies 
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FIG. 1. A nonsymmorphic symmetry {g|t} leads to band 
crossings on a g invariant line or plane in momentum space, 
(a) Without other symmetries, pairs of bands intersect an odd 
number of times as they cross the BZ. (b) With time-reversal 
symmetry, with 0^ = -|-1, the crossing occurs at the zone 
boundary G/2, where = —1. For 0^ = —1, Kramers 

pairs at k = 0 and G/2 connect as in (c), leading to a line 
node in an invariant plane or a Weyl node on an invariant line. 
The labels indicate the eigenvalues of (d) With 

inversion and time-reversal symmetry all states are degenerate 
(they are offset for clarity), and the crossing occurs at G/2. 


pairs of band branches must intersect an odd number of 
times as they cross the BZ, as shown in Fig. 1(a) . 

Time-reversal symmetry 0 imposes further con¬ 
straints. The situation is simplest in the absence of spin- 
orbit interactions, where effectively 0^ = -1-1, and = 1, 
so A = ±1. In this case the degeneracy must occur at the 
time-reversal invariant momentum k = G/2, at the BZ 
boundary [Fig 1(b)]. At that point the eigenstates of 
{^jt} with eigenvalue ±i are interchanged by 0. Stated 
another way, at k = G/2, the Bloch Hamiltonian com¬ 
mutes with 0 = {(;|t}0, which satisfies 0^ = —1, guar¬ 
anteeing a Kramers degeneracy. 

Spin-orbit interactions lead to additional splitting of 
the bands, though time-reversal symmetry (with 0^ = 
— 1) enforces Kramers degeneracies at the time-reversal 
invariant momenta. Moreover, a mirror or twofold rota¬ 
tion satisfies = —1, so that {gjt} has eigenvalues 
at k = 0, and eigenvalues ±1 at k = G/2. It follows that 
Kramers partners at k = 0 have opposite eigenvalues un¬ 
der {gjt}, while Kramers partners at k = ±G/2 have 
the same eigenvalue under {gjt}. In the absence of other 
symmetries, this leads to the pattern of splitting shown 
in Fig. 1(c). This results in the sticking of four bands, 
so that a system with the nonsymmorphic symmetry and 
time-reversal that has a band filling of 4n -I- 2 (for integer 


n) is necessarily a topological semimetal. Note that for 
a glide mirror plane g = this leads to a line of degen¬ 
eracies, while for a glide mirror line g = Mn^ , it leads to 
(Weyl) point degeneracies. 

If in addition the crystal has inversion symmetry P, 
then since (P0)^ = —1, the bands are Kramers degener¬ 
ate for all k. This leads to a fourfold degenerate crossing 
at k = G/2 [Fig. 1(d)]. At k = G/2, P and {(/[t} both 
have eigenvalues ±1 and commute with 0. In addition 
they anticommute with each other, which guarantees a 
fourfold degeneracy. Since the fourfold degeneracy will be 
split away from G/2, this constitutes a 2D Dirac point. 

We now introduce a simple tight-binding model for 
a 2D Dirac semimetal. This can be viewed as a 2D 
analog of the diamond lattice model for a 3D Dirac 
semimetal [21]. We begin with a square lattice of s states 
with first and second neighbor hopping and create a 
X unit cell [Fig. 2(a)]. This features line nodes 
along the BZ boundary. To make the lattice nonsym¬ 
morphic we displace the A{B) sublattices in the +z{—z) 
direction. This crinkling of the lattice [Fig. 2(b)] permits 
a second neighbor spin-orbit interaction [8], leading to the 
Hamiltonian 

k k 

H =2tTx cos cos + t 2 (cos kx + cos ky) 

2 2 

Pt^'^TzlaySinkx-CTxSvnky], ( 1 ) 

where r and a are Pauli matrices describing the lattice 
and spin degrees of freedom. As shown in Fig. 2(b), 
nonzero breaks the degeneracy on the zone boundary 
everywhere except the corners M and edge midpoints Xi 
and X 2 , at which appear Dirac points. 

The present lattice has high symmetry (layer group 
PA/nmm) with multiple symmetries protecting the Dirac 
points. From the analysis in Fig. 1(d), the combination of 
P, 0, and the glide mirror plane symmetry {Mz\W} (t is 
in units of the Bravais lattice constant) protects the Dirac 
points at Xi and AI 2 , while the screw axes {C 2 x\^ 0 } and 
{C 2 y\Q\} protect the Dirac points at Xi,M and X 2 ,M, 
respectively. This can be further seen by examining the 
k ■ p theory near these points. Near k = M, 

niM + q) = t^'^iT^ayQx - T^Gxqy). ( 2 ) 


At M, the symmetries 0 = iayK and P = Tx allow 
a single mass term Tx- This is forbidden by {C 2 x\h^} = 
Tyax and {C 2 y\Q\} = Tyay, but is allowed by = 

iTx<yz- Likewise, near k = Xi, 


'H(Xi +q) = {tTx - P T^(Ty)qx + t 




cqy 


( 3 ) 


At Xi, 0 = iTzTyK and P = Ty allows the mass terms Ty. 
This is forbidden by {^241^0} = Txax and = 

TxO'z, but is allowed by {C2y\0^} = iTyay. A similar 
analysis applies to X 2 . 

In addition to spatial symmetries, Eq. (1) exhibits 
a particle-hole symmetry when t 2 = 0, manifested by 
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FIG. 2. Energy bands for structures with Dirac points pro¬ 
tected by inversion symmetry, with lattice structure and BZ 
shown above. Dirac points and nodal lines in the BZ are 
marked in cyan, (a) The x supercell of a square lattice 
with nodal lines along the BZ edge protected by (b) 

Crinkling the lattice breaks {i5|t}, leaving nonsymmorphic 
symmetries and Dirac points at Xi, X2, and M. (c) Distort¬ 
ing in the (11) directions eliminates screw axes {C 2 a;|| 0 } and 
{C2y\0^} (as well as C4), gapping the corner Dirac point, (d) 
Alternatively, displacing one of the sites along (10) breaks 
leaving Dirac points at the corner and one edge. 


{7i,Ty} = 0. This guarantees the Dirac points at Xi, 
X 2 , and M occur at the same energy. t 2 violates this 
symmetry and leads to a shift in the energy at M rel¬ 
ative to Xi, 2 - However, the mirror lines Mx±y guaran¬ 
tee the equivalence of Xi 2 . Thus, while the touching 
of the conduction and valence band at Xi, X 2 , and M 
is guaranteed, there will in general be electron and hole 
pockets with a finite Fermi surface [Figs. 2(b) and 2(d)]. 
Nonetheless, with appropriate band structure engineer¬ 
ing it may be possible to tune the edge and corner Dirac 
points to approximately the same energy. In the fol¬ 
lowing we will explore the range of behaviors that can 
arise when the symmetries in Eq. (1) are systematically 
lowered. We show that Eq. (1) lies at the boundary of 
three distinct Dirac semimetal phases, two with a pair of 
Dirac points and the third with three, and prove that in 
2D there cannot exist a single symmetry-protected Dirac 
point. 

Case I: Two symmetry equivalent Dirac points- 
Eirst, we consider a distortion that breaks the sym¬ 
metry between interactions in the (11) and (11) direc¬ 
tions [Eig. 2(c)], but preserves the mirror line M^+y. This 


allows a distortion of the first neighbor hopping term, 

Vi = Ai sin y sin (4) 

The Hamiltonian H + Vi retains inversion P, along with 
and {M^j+yjOO}, but the C 2 screw symmetries 
are violated. As shown in Fig. 2(c), the corner Dirac 
point is gapped, but the Dirac points at 2 remain. 
Provided there are no other extraneous electron and hole 
pockets and the electron count is 4n -b 2, a system with 
these symmetries will be a 2D Dirac semimetal, with two 
symmetry equivalent Dirac points at Ep- In fact, this is 
the only truly protected Dirac semimetal in the absence 
of approximate particle-hole symmetry, ^2 0. If M^+y 

is violated, then the 2 Dirac points are inequivalent, 
which belongs in the next case. 

Case II: Two symmetry inequivalent Dirac points.-We 
next consider breaking the glide mirror plane 
while keeping {C'2il|0}. We displace the B site in the y 
direction [Fig. 2(d)], allowing a term 

1/2 = A2COSy sin^Ty. (5) 

The Hamiltonian H + V 2 now has a gap at X 2 , but Dirac 
points remain protected at Xi and M, though in the 
presence of t 2 they are at different energies. It follows 
from the arguments illustrated in Fig. 2 that a glide plane 
with fractional translation (^O) will produce this result as 
well. 

Case III: Three Dirac points.-AW three Dirac points 
remain protected in the presence of {C' 2 £l^ 0 }, 

and {C2y\0^}- We find that the three Dirac points can 
persist even when the screw symmetries are violated pro¬ 
vided the system retains an additional Ca symmetry 
about the center of a plaquette. This is violated by the 
crinkling responsible for t®®. However, for the flat sys¬ 
tem tso = 0 this symmetry will pertain if the A and B 
sites have a lower symmetry internal structure allowing 
a spin-orbit term 

^3 = t'soDi- sin k^aj: + sin kyay) (6) 

that preserves the three Dirac points. Note, however, 
that it is impossible to gap the Xi and X 2 points without 
also gapping M. As discussed below, a single symmetry 
protected Dirac point is not possible. 

These three cases represent all possible 2D systems 
with Dirac points: (I) symmetric Dirac points at both 
edges protected by P and a glide plane, (2) Dirac 
points at an edge and corner protected by P and glide 
planes or screw axes along the x or y directions, and (3) 
Dirac points at both edges and the corner protected by 
P and a C^z rotation, or case H for both the x and y 
directions. 

Case IV: Line nodes and Weyl pomts.-Finally, we 
mention that if inversion symmetry is violated, while 
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FIG. 3. Without inversion, nonsymmorphic symmetries pro¬ 
tect nodal lines or Weyl points. This may be achieved by 
altering each site to be bipartite and asymmetric (e.g., het¬ 
erodimers). (a) Breaking inversion in case I, while preserving 
leads to nodal lines that circle the former Dirac 
points at Xi^ 2 - (b) Breaking inversion in case II, preserving 
{C2x\^0}, results in Weyl points on the (724 invariant lines. 


keeping the nonsymmorphic symmetries, then the Dirac 
points are removed, but there remain Weyl points or line 
nodes. For example, if in case I each site has a dipole mo¬ 
ment p(x + y), then a spin-orbit term vsqTxO'z si'ii{kx — 
ky)/2 is allowed. This preserves and for small 

uso leads to a circular line node surrounding the erst¬ 
while Dirac points at Xi ^2 [Fig. 3(a)]. More generally for 
this symmetry, a line node will separate regions contain¬ 
ing F and M from regions containing Xi^ 2 - Similarly, if 
{C 2 xW^} is preserved but not then the two 

Dirac points each split into two Weyl points along the 
C 2 x invariant lines ky = 0 and ky = tt [Fig. 3(b)]. Weyl 
points at generic k are also locally protected when 0 
and a C 2 z symmetry are preserved even in the absence of 
nonsymmorphic symmetries. However, in this case the 
Weyl points can annihilate, and are not guaranteed by 
symmetry. 

A distinctive feature of the nonsymmorphic Dirac 
semimetals is that they describe a critical point separat¬ 
ing topologically distinct phases. By lowering the sym¬ 
metry it is possible to open a gap that leads to either a 
trivial or topological 2D insulator. Consider the system 
of case I [Fig. 2(c)] with a gap introduced by displacing 
the center atom [Fig. 4(a)], described by a perturbation 


V4 


mi sin 



+ m 2 sin 



where mi and m 2 describe the displacement in the [11] 
and [11] directions, respectively. When jmij > |m 2 |, the 
system is a topological insulator, and for jmij < |m 2 | 
the system is in a trivial phase [Fig. 4(b)]. The bound¬ 
aries between topological and trivial insulating phases 
are marked by a nonsymmetry-protected Dirac point. 
Changing the sign of mi or m 2 results in a phase with 
the same topological character. However, one may gap 
the Dirac semimetal directly into either a topological or 
trivial insulator, depending on the direction of the dis¬ 



(a) 



FIG. 4. By further lowering the symmetry the Dirac 
semimetal in Fig. 2(c) can be driven into either trivial (I) 
or topological (TI) insulator phases. Lattice displacements 
in (a) give rise to the phase diagram in (b). 


placement. 

It is worthwhile to contrast the 2D Dirac semimetal 
presented here with the symmetry protected Dirac 
points that arise at the surface of a weak topological 
insulator[21] or a topological crystalline insulator[22]. 
There, the surface Dirac points are protected by the com¬ 
bination of time-reversal and a translation or mirror sym¬ 
metry. Breaking the symmetry leads to topologically dis¬ 
tinct gapped surface phases that map to each other under 
the symmetry. Importantly, this leads to an absence of 
localization when the surface is strongly disordered, but 
retains the symmetry on the average[23-26]. This phe¬ 
nomenon can occur only at the surface of a 3D topological 
phase. In a purely 2D system the trivial and TI phases 
are not related by symmetry, so our 2D Dirac semimetal 
can be localized by strong disorder. However, since it is in 
the symplectic class, weak disorder leads to antilocaliza¬ 
tion. The absence of symmetry relating the trivial and 
TI phases rules out a single symmetry protected Dirac 
point, since in that case the symmetry would change 
the sign of the single mass term. We find that chang¬ 
ing the sign of the mass of one of our Dirac points by 
a symmetry breaking perturbation always leads to the 
change in sign of another Dirac point, resulting in the 
same topological order. Weak electron-electron interac¬ 
tions do not significantly alter the electronic structure of 
a 2D Dirac semimetal, though strong interactions could 
lead to a gapped state. However, if the symmetry is not 
lowered that state must exhibit a nontrivial topological 
order [20]. 

In terms of realistic materials, there are both 
challenges and advantages to working in 2D. Two- 
dimensional systems can be more fragile, and typically 
require substrates that can influence the behavior. On 
the other hand, they offer additional tunability not avail¬ 
able in 3D systems. For example, the proximity of TI and 
I phases could allow for flexible patterning of helical edge 
channels. There has been surge of interest in monolayers 
of covalently bound layered materials over the past few 
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years [27]: these materials exhibit abundant variety in 
composition and structure [28], and have already shown 
hints of being able to host graphenelike Dirac points [29] 
and topological phases [30, 31]. Many of these materi¬ 
als exist in structures belonging to appropriate symme¬ 
try groups, including litharge (case III) and WTe 2 [32] 
(case II), and some are known to possess the appropri¬ 
ate band filling, such as (Nb,Ta)Te 2 [33]. Moreover, the 
Dirac points in iridium oxide superlattices proposed by 
Chen and Kee[18] constitute a manifestation of case I de¬ 
scribed above. They showed that distortions can lead to 
a TI gapped phase. It will also be interesting to demon¬ 
strate the distortions in that system that lead to the triv¬ 
ial insulator. We are thus optimistic about the prospects 
for the experimental study of 2D Dirac semimetals. 
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